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ABSTRACT 

We present 0.15" resolution observations of the 227 GHz continuum emission from the circumstellar 
disk around the FU-Orionis star PP 13S*. The data were obtained with the Combined Array for 
Research in Millimeter- wave Astronomy (CARMA) Paired Antenna Calibration System (C-PACS), 
which measures and corrects the atmospheric delay fluctuations on the longest baselines of the array 
in order to improve the sensitivity and angular resolution of the observations. A description of the 
C-PACS technique and the data reduction procedures are presented. C-PACS was applied to CARMA 
observations of PP 13S*, which led to a factor of 1.6 increase in the observed peak flux of the source, 
a 36% reduction in the noise of the image, and a 52% decrease in the measured size of the source 
major axis. The calibrated complex visibilities were fitted with a theoretical disk model to constrain 
the disk surface density. The total disk mass from the best fit model corresponds to 0.06 M Q , which is 
larger than the median mass of a disk around a classical T Tauri star. The disk is optically thick at a 
wavelength of 1.3 mm for orbital radii less than 48 AU. At larger radii, the inferred surface density of 
the PP 13S* disk is an order of magnitude lower than that needed to develop a gravitational instability. 
Subject headings: stars: circumstellar matter — stars: individual(PP 13S*) — stars:pre-main sequence 
— techniques: interferometric 



1. INTRODUCTION 

Electromagnetic waves from an astronomical radio 
source suffer distortion f rom irregularities in th e refrac- 
tivity of the atmosphere ( Thompson et al.|2001 ). At mil- 
limeter and submillimeter wavelengths, the distortions 
are caused primarily by a turbulent water vapor distri- 
bution, though dry air tur bulence may also be important 
under some circumstances (Stirling et al.|2006). The spa- 



the effect of atmospheric delay fluctuations 
forward 



tial scales of the turbulence extend to kilometer distances 
with a power law spectrum ( |Kolmogorov| 1 94l| ) that cre- 
ates atmospheric delay fluctuations on timescales that 
range from fractions of a second to tens of minutes. The 
signal degradation is particularly serious for millimeter- 
wave radio interferometers in extended configurations, 
where perturbation of phases across the instrument of- 
ten exceeds many radians. These perturbations can lead 
to decorrelation (loss of am plitude), d istortion of the im- 
age, and loss of resolution (Lay 1997). 
Different approaches have been employed to overcome 

1 Department of Astronomy, California Institute of Technol- 
ogy, 1200 East California Blvd, Pasadena, CA 91125 

2 Owens Valley Radio Observatory, California Institute of 
Technology, Big Pine, CA 93513 

3 Department of Astronomy, University of Maryland, College 
Park, MD 20742-2421 

4 Combined Array for Research in Millimeter-wave Astron- 
omy, P.O. Box 968, Big Pine, CA 93513 

5 Department of Astronomy and Astrophysics, University of 
Chicago, 5640 S. Ellis Ave. Chicago, IL 60637 

6 Space Sciences - VP62, NASA Marshall Space Flight Center, 
Huntsville, AL 35812 

7 Department of Astronomy, University of Illinois, Urbana, IL 
61801 

8 Jansky Fellow, National Radio Astronomy Observatory 

9 Astronomy Department, University of California at Berkeley, 
Berkeley, CA 94720-3411 



A s t raight - 
approach is self- calibration (Schwab 1980| , 
where the visibility phase used to calibrate the data is 
measured from the actual science target. A model for 
the source spatial structure is needed and bright sources 
are required to measure the fringe phase with a high sig- 
nal to noise ratio in a short integration. Two other ap- 
proaches have bee n applied to faint targets. In water 



vapor radiometry ( Westwater|| 1967 Woody et al.||2000 ) 
a dedicated radiometer monitors the water vapor emis 



le water vapor emis- 
sion along the poin ting direction of the a ntenna. In fast 



position switching (Holdaway et al. 1995), the antennas 



switch rapidly between the science target and a nearby 
phase calibrator to capture the atmospheric delay fluc- 
tuations on time scales longer than the switching cycle 
time. Both methods probe the atmosphere close to the 
line-of-sight toward the science target, but have limita- 
tions in actual applications. Fast position switching re- 
duces the time spent on source by a factor of ~ 2, while 
water vapor radiometry requires a physical model to re- 
late the water line brightness and the path correction, 
as well as it assumes that refractivity fluctuations are 
produced only by water vapor. 
As an al ternative approach , an a rray of closely paired 



antennas (Asaki et al. 19961 1998) continuously moni 



tors the atmospheric phase fluctuations by observing a 
nearby calibrator. Two arrays of antennas are neces- 
sary: antennas belonging to the "science array" observe 
the science target and phase calibrator, while antennas 
belonging to the "reference array" simultaneously mon- 
itor an atmospheric calibrator. Phase correction on the 
science target and phase calibrator is accomplished by 
subtracting the visibility phase measured from the at- 
mospheric calibrator on each baseline. An advantage of 
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the paired antenna technique over water vapor radiome- 
try is that it accounts for atmospheric phase fluctuations 
due to both water vapor and a dry air component. 

Between November 2008 and February 2009, the Com- 
bined Array for Research in Millimeter-wave Astronomy 
(CARMA) implemented the CARMA Paired Antenna 
Calibration System (C-PACS) to correct for atmospheric 
delay fluctuations on the longest baselines of the inter- 
ferometer (up to 2 km in length). The goal of C-PACS is 
to enable routine imaging in the most extended CARMA 
configurations. 

In this paper we describe C-PACS and apply this cal- 
ibration technique to observations of the circumstellar 
dust around PP 13S =t j^2l PP 13S* is a young pre-main 
sequence star located in the constellation of Perseus and 
em bedded in the L 1473 dark cloud at a distance of ~ 350 
pc ( |Cohen et al. [1983 ). This object is thought to have 
experienced a FU-Orionis- type outburst in the past due 
to a massive accretion episode and is now declining in 



brightness to a quiescent state (Aspin & Sandell 
The FU-Orionis nature of PP 13S* has been esta 



2001J). 
}lished 



based on the highly broaden ed infrared CO absorption 
bands ( |Sandell fc Aspin|1998[ ), the jet-like feature seen in 
[SII] e mission which is characteristic of Herbig-Haro out - 
flows ( |Aspin fe Reipurth|[2050) JAspin fc Sandell||2001| ), 



and the consistent dimming and morphology changes ob 
served a t near-IR and optical w avelengths over several 
decades ( |Aspin fc Sandell||2001[ ) . All of th ese character- 
istics are com mon to FU-Orionis objects ( jHartmann fc| 
Kenyon | 1996|. The star, with a bolometric luminosity of 
30 L ( Cohen et al.|1983 ), is surrounded by an extended 
disk and envelope that conta ins about 0.6 M of gas and 
dust ( |Sandell fc Aspin||1998[ ). The new CARMA obser- 
vations will help understand the origin of the FU-Orionis 
phenomena in PP 13S* by measuring the distribution of 
dust continuum emission on spatial scales of ~ 50 AU. 

2. DESCRIPTION OF C-PACS 

Before presenting the new observations of PP 13S*, we 
describe the paired antenna calibration system as imple- 
mented at CARMA. We first describe the antenna con- 
figuration used for the observations, and outline the basic 
principles of the technique. 

2.1. CARMA 

CARMA is an heterogeneous interferometer compris- 
ing 23 antennas: six 10.4-m telescopes from the Califor- 
nia Institute of Technology/ Owens Valley Radio Obser- 
vatory (OVRO), nine 6.1 m telescopes from the Berkeley- 
Illinois-Maryland Association (BIMA), and eight 3.5- 
m telescopes from the University of Chicago Sunyaev- 
Zel'dovich Array (SZA). The two most extended array 
configurations contain baselines that range in length from 
100 m to 1000 m (B configuration) and 250 m to 1900 m 
(A configuration) to achieve an angular resolution of 
0.3" and 0.15" respectively, at an observing frequency 
of 230 GHz. 



A schematic of C-PACS is shown in Figure [T] C-PACS 
pairs the eight 3.5-m antennas with selected 6-m and 10- 
m antennas, preferentially on the longer baselines. Typ- 
ically, a 3.5-m antenna is offset by 20-30 m to the west 
of a larger antenna. This separation is a compromise 
between the need to put the antennas as close as pos- 
sible to probe the same atmospheric path and to avoid 
shadowing between antennas. 

The science array, composed of the 6-m and 10-m an- 
tennas, operates in the 1.3- or 3-mm atmospheric win- 
dows as requested by the investigator. The reference 
array, comprising the 3.5-m antennas, operates in the 
1-cm window. The observing cycle consists of observa- 
tions of the science target interleaved with periodic ob- 
servations of the phase calibrator. Both the science and 
reference arrays observe the phase calibrator to measure 
instrumental phases drifts. Subsequently, the science ar- 
ray observes the science target while the reference array 
monitors a strong point source (i.e., the "atmospheric" 
calibrator) close to the science target, to measure the 
delay introduced by the atmosphere. The atmospheric 
delay measured by the 3.5-m antennas can then be ap- 
plied to the science observations. 

2.2. Properties of the atmosphere 

It is helpful to have a physical picture to understand 
the principles of the correction. We suppose the atmo- 
sphere to be a pattern of random refractive index varia- 
tions that is blown across the array at the wind velocity. 
Furthermore, we assume that the layer is at a height of a 
couple of kilometers and that the thickness is much less 
than the height. General experience at this and other 
sites show that these conditions are often consistent with 
what is observed. Other conditions can be present, but 
can often be characterized by two or more layers at dif- 
ferent altitudes with separate wind vectors so that only 
a small change to the analysis is required. 

Kolmogorov theory predicts a turbulence distribution 
with a power law spectrum from less than a millimeter 
in size to many kilometers. This results in random delay 
differences between signals arriving at different antennas 
that increase w ith separation as a power law function 
(Tatarskii 1961). As the pattern moves over the array, 
the delay differences are observed as temporal fluctua- 
tions in the visibility phases. The RMS of the delay 
depends on the wind speed, but not its direction. Struc- 
tures smaller in size than the antenna diameter are av- 
eraged out and do not contribute to the phase errors. 
Structures on scales large in comparison with the base- 
line length are common to the two antennas on the base- 
line and therefore tend to cancel out. From these the- 
oretica l considerations sup ported by experimental evi- 
dence ( |Sramek|1983 19891, it is found that the resulting 
delay variance vs baseline length (i.e., the delay struc- 
ture function) also follows a power-law. The theoretical 
slope of the power-law is 5/3 and 2/3 for two- and three- 
dimensional Kolmogorov turbulence, respectively. 



10 pp 23 a come tary nebula in the list of |Parsamian fe Pet-| 
|rossian| ( |1979| l. PP 13S is the southern component containing a red 
ne bula with a bright in frared point source at the apex as designated 
bylCohen et al.|(|1983k. PP 13S* corresponds to the embedded star 
itself, which is obscured by circumstellar material. The northern 
component, PP 13N, is a T Tauri star. 



2.3. Atmospheric delay corrections 



Following the discussion in Asaki et al. ( 1996 ) , consider 
two pairs of antennas as shown in Figure [I] In the sim- 
plest case, we measure the delay difference between the 
antennas on the reference baseline (i.e., baseline A3 — A± 
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in Figure [I]) as a function of time. Assuming a non- 
dispersive atmosphere, the delay on the science baseline 
(i.e., baseline A\ — A 2 ) is corrected by applying the de- 
lay difference on the reference baseline to the visibility 
measurement. 

The reference and science delays are not identical since 
the two baselines are not exactly co-located. The relevant 
distances that determine the efficacy of the delay correc- 
tions are not the baseline lengths at the ground, but the 
distances between the radio beams as they traverse the 
turbulent layer (i.e. A' x — A' 3 and A' 2 — A' 4 in Figure 
[T]). The beam separation at the turbulent layer depends 
upon the relative positions of the target and reference 
source in the sky, the height of the turbulent layer, and 
the configuration of the antennas on the ground. The 
upper limit to the beam separation is given by 



the measured phase from radiometer noise is given by 



\Ai — A2I + ah/sin(e), 



(1) 



A' 4 will be much less than either 
A' 2 or the reference beam 



where a is the angular separation between the science 
target and the atmospheric calibrator, h the height of 
the turbulence, and e is the source elevation. Assuming 
that the turbulent layer is at a height of 1 km (contin- 
uous line) or at 2 km (dashed line), Figure [2 shows the 
trajectories of the 3.5-m beam locations relative to the 
6-m and 10-m beams for an 8 h observation of PP 13S* 
centered on transit with 3C111 as the atmospheric cali- 
brator. As shown in this figure, the choice of 3C111 as 
the atmospheric calibrator is particularly fortuitous for 
PP 13S* since the science and reference beams for most 
antennas nearly cross wi thin the turbule nt zone. 

Using phase closure (Jennison 1958), the difference 
(At) between the actual delay for the target beams 
A\ — A' 2 and the atmospheric reference beams A' 3 — A' 4 is 
given by At = t(A[ - A' 2 ) - t(A' 3 - A' A ) = r(A[ - A' 3 ) - 
t(A' 2 — A' 4 ). In a favorable configuration, the beam sepa- 
rations A[ — A' 3 and j 
the target beam separation A\ 

separation A' 3 — A 4 . The RMS of the corrected visibilities 
will be V2 worse than for an array with beam separation 
A\ — A 3 , which implies that C-PACS will have a perfor- 
mance equivalent to an array that has baseline lengths 
~20-70% larger than the A\ — A 3 beam separation, de- 
pending on the structure function exponent. The com- 
plete analysis contains additional correlation terms and 
added uncertainties caused by the finite signal-to-noise 
for the atmospheric calibrator observations. 

2.4. Atmospheric calibrators 

The ability of C-PACS to correct the atmospheric de- 
lays is limited by the delays from the short beam spacings 
A[ — A' 3 and A' 2 — A 4 (see Figure I]) , instrumental phase 
drifts on the reference array, andthe radiometer noise. 
The delay errors caused by differences in the beam spac- 
ings between the science and reference arrays are given 
statistically by the structure function R{\A\ — A'j\). The 
instrumental errors can be removed by removing a box- 
car average over the length of the observation, and will 
contribute negligible delay errors as long as the timescale 
for the instrumental drifts are large compared to the box- 
car width. The delay errors due to radiometric phase 
noise depend on the strength of the source being ob- 
served, the receiver properties, and the atmospheric char- 
acteristics (Thompson et al. 20011. The uncertainty in 



Acj) 
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where fee is the Boltzmann's constant, T sys is the system 
temperature, r/Q is the correlator quantization efficiency, 
A e s is the effective collecting area of the antennas, S is 
the flux density of the atmospheric calibrator, B is the 
bandwidth of the observations, and t is the integration 
time. The measurement uncertainty in the delay is then 
A(f>/ \2irv) . The net variance in the delay after applying 
C-PACS is given by 



A7 



k B T, 



B-Lsys 



^■q Q A c ttSiiv\/2BtJ 

The structure function is generally described by a power 
law with exponents varying from 5/3 to 2/3 depending 
upon the spacing and thickness of the turbulent layer. 
The scaling coefficient of the power law also varies de- 
pending upon the weather conditions. In order for C- 
PACS to improve the image quality, the target and ref- 
erence beams must be close at the turbulent layer such 
that R(\A[ - A 3 \) + R(\A 2 - A 4 \) < A 2 . This 'requires 
angular separations < 5° for the A and B configuration 
C-PACS pairings and typical winter weather conditions 
at the CARMA site. A future publication will use ac- 
tual measurements to quantify how the quality of the 
C-PACS correction varies with angular separation be- 
tween the science target and the atmospheric calibrator 
(Zauderer et al., in preparation). 

The radiometer noise should also contribute much less 
than a wavelength of delay error for the C-PACS cor- 
rections to be successful. For the characteristics of the 
3.5-m telescopes and the 1-cm receivers, the radiometer 
delay error is given by 



AT ra dj, 
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VUy/ 



Vis 



-1/2 



(4) 



Thus, 1.3-mm observations with integration times of 
t = 4 s (short enough to measure and correct most of 
the atmosphere fluctuations) require a reference source 
brighter than S ~ 1 Jy in the 1 cm band. When sev- 
eral atmospheric calibrators are available, the optimum 
choice between calibrator separation and brightness can 
be found by minimizing Equation [3] for the expected 
weather conditions. 

We combined the SZA 30 GHz calibrator list the GBT 
catalog at 1.4 and 5.0 GHz dCondon fc Yin 200 l|, a nd the 
WMAP point source catalog (f\\ right et al.||2009p to es- 
timate the density of potential C-PACS calibrators. For 
each source in the GBT catalog, we extrapolated the flux 
density from 5.0 GHz to 30 GHz by measuring the spec- 
tral index a between 1.4 and 5.0 GHz (S v cx v a ). We find 
that 50% of the sky is within 5° of a point source with 
flux density greater than 1 Jy at 30 GHz. The num- 
ber of suitable C-PACS calibrators could be expanded 
by increasing the sensitivity of the reference array, which 
would allow us to employ fainter atmospheric calibrators. 
This could be accomplished by increasing the correlator 
bandwidth or improving the receiver sensitivity. 
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3. OBSERVATIONS AND DATA REDUCTION 

PP 13S* is particularly well suited for C-PACS obser- 
vations since the nearest atmospheric calibrator (3C111) 
is bright (~4 Jy at 1.3 mm at the time of the obser- 
vations) and separated by 1.5° from PP 13S*. Thus the 
calibrator satisfies the basic criteria needed for successful 



C-PACS corrections (see Section 2.4). In this section, we 



describe the CARMA observations and data reduction of 
PP 13S*. 

3.1. CARMA 1.3-mm wavelength observations 

The 6-m and 10-m antennas were used to obtain 
1.3 mm continuum observations of PP 13S* on UT Dec 5, 

2008 in the CARMA B configuration and on UT Jan 18, 

2009 in the CARMA A configuration. Double-sideband 
receivers mounted on each antenna were tuned to a rest 
frequency of 227 GHz placed in the upper sideband. The 
correlator was configured with three 468.75 MHz wide 
bands to provide 1.41 GHz of continuum bandwidth per 
sideband. The observing sequence interleaved 3 min ob- 
servations of 3C111 with 12 min observations of PP 13S* 
in B configuration and 4 min in A configuration. The 
complex visibilities were recorded every 4 s. 

Data reduction was performed using the Multichan- 
nel Image Reconstruction, Imag e Analysis and Display 



(MIRIAD) software package (Sault et al. 



19951. Each 



night of observations was calibrated separately. The cal- 
ibration consisted of first applying a line-length correc- 
tiorpl and then a passband correction derived from ob- 
servations of 3C111. Only these two calibrations were 
applied to the millimeter data befor e pro ceeding with 
the C-PACS corrections (see Section 3.3). All images 



that are presented here were formed by inverting the vis- 
ibility data using natural weighting, and then "cleaning" 
with the point spread function (i.e. the "dirt y" beam) 
using a hybrid Hogbom/Clark/Steer algorithm ( jHogbom] 
1974] |Clark||1980l |Steer, Dewdney, fc Ito||1984[ ). 



3.2. CARMA 1-cm wavelength observations 

The eight 3.5-m antennas were used to obtain 1-cm 
observations of the atmospheric calibrator simultane- 
ously with the 1.3-mm wavelength observations. Single- 
sideband receivers mounted on each antenna were tuned 
to a sky center frequency of 30.4 GHz. For these an- 
tennas, a wideband correlator is available that was con- 
figured with fourteen 500 MHz wide bands to provide 
7 GHz of continuum bandwidth. Complex visibilities 
were recorded every 4 s in order to track the rapidly 
varying atmospheric fluctuations. Single-sideband sys- 
tem temperatures for the 1-cm observations ranged be- 
tween 35 and 55 K. 

The data calibration consisted of applying a time- 
dependent passband measured from an electronically cor- 
related noise source that was observed every 60 s. This 
passband was computed and applied on a 60 s timescale 
to remove any delay variations in the digitizers due to 
temperature cycling of the air-conditioning that cools the 
correlator. A passband correction derived from observa- 
tions of 3C111 was then applied. 

11 The CARMA line-length system measures the total round 
trip delay caused by possible mechanical effects and temperature 
variations of the fiber-optic cables running to each antenna. 



3.3. Applying 1-cm delays to the 1.3-mm data 

Phase referencing would normally be performed at this 
point in the data reduction process to remove the slowly 
varying phase drift introduced by the instrument. How- 
ever, a phase calibration computed over a long time inter- 
val, and prior to correcting for the "fast" atmospheric de- 
lay fluctuations will alias the fast compo nent into a sl owly 
varying error on the phase calibration (| Lay||1997b ). C- 
PACS can reduce the errors introduced from standard 
phase calibration techniques by correcting for both fast 
and slow atmospheric delay fluctuations. 

The delays were extracted from the 1-cm wavelength 
observations of the atmospheric calibrator, and then ap- 
plied to each corresponding paired baseline in the science 
array. Only eight of the 15 antennas from the science 
array, those paired with a reference antenna, have the 
C-PACS correction applied. The delay derived from the 
reference array is applied to the science target (PP 13S*) 
and phase calibrator (3C111) on each record. The delays 
were computed from the mean phase across all channels 
in the 1-cm observations divided by the mean frequency. 
The observed wavelength of 1 cm used for the reference 
array is longer than the typical atmospheric delay fluctu- 
ations, and delay tracking through phase wraps was not 
a problem. 

The atmospheric delays derived from the 1-cm data 
were applied directly to the 1.3-mm data without any 
corrections for differences in the observed frequency. 
This is possible since the dispersion in refractivity of wa- 
ter vapor between centim eter and m illimeter wavelengths 
is less than a few percent ( |Hill|1988[ ) away from the strong 
atmospheric emission lines, an d ionospheric effec ts are 
negligible at these frequencies ( |Hales et al.| 2003) . The 
minimal dispersion of the refractivity has also been ver- 
ified ex perim entally from the C-PACS observations (see 
Section 4.1). After applying the delay corrections to the 
1.3 mm data on 4 s intervals, a long-interval (10 min) 
phase calibration is applied to the 1.3-mm data to re- 
move the slow varying instrumental delay difference from 
the two different arrays. 

4. APPLICATION OF C-PACS 

The effect of the C-PACS on the calibrated phases can 
be analyzed in several stages. We first compare the phase 
fluctuations measured toward 3C111 at wavelengths of 
1 cm and 1.3 mm to demonstrate that the reference and 
science arrays are tracing the same atmospheric fluctua- 
tions. We then demonstrate that C-PACS yields quan- 
titative improvement in the quality of the PP 13S* star 
image. Throughout this section, no absolute flux scale or 
amplitude calibration are applied to the data in order to 
evaluate how C-PACS im prov es the phase stability o n th e 
phase calibrator (Section 4T]) and PP 13S* (Section |42|). 



4.1. 3C111 

Figure [3] shows the measured fringe phase toward 
3C111 for one paired baseline in the B configuration to 
illustrate the correlation that exists between the phases 
measured at wavelengths of 1.3 mm and 1 cm. For this 
figure, the phases measured on the reference array were 
scaled by the ratio of the observed frequencies (227 GHz 
/ 30.4 GHz = 7.5). The correlation between the 1-cm 
and 1.3-mm phases is evident over the nearly 8 h time 
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period and is present for all paired baselines. In ad- 
dition, the phases are tracked between the two arrays 
even though the science array is switching between two 
sources. These results demonstrate that (1) the 1-cm 
phases can be used to track the delay fluctuations at 
higher frequencies, (2) the atmosphere is non-dispersive 
at these wavelengths such that a linear scale with fre- 
quency can be used to predict the phases fluctuations at 
other wavelengths, and (3) C-PACS can correct the sci- 
ence observations while preserving the link to the phase 
calibrator observations in the science array. 

The observed phases at 227 GHz toward 3C111 on all 
paired baselines in the A and B configurations are shown 
as a function of uu-distancf ^1 in Figure H The upper 
panels show the visibility phases measured in 4 s inte- 
grations before applying C-PACS corrections, and the 
middle panels show the visibility phases after applying 
the corrections. The RMS scatter on the visibility phase 
before correction is between 33° (at short itv-distances) 
and 53° (at long ww-distances) for B configuration and 
between 26° and 48° for A configuration. After apply- 
ing the C-PACS correction, the phase scatter is reduced 
to 15 — 18° for correction across all wu-distances in both 
configurations. 

Another way to grasp the effect of the C-PACS cor- 
rections can be seen on the bottom panels in Figure [4j 
where the coherence value (e~* RM s/ 2 ) is measured over 
uw-distance bins of width 130 kA. Before applying the 
C-PACS corrections, the coherence decreases with in- 
creasing baseline length since the longer baselines have 
larger atmospheric fluctuations. After applying the C- 
PACS corrections, the coherence is uniform with base- 
line length at a value of ~ 95%. The coherence becomes 
nearly constant with baseline length since C-PACS con- 
verts the 200-1800 m baselines into ~ 30-50 m effective 
baselines for all paired antennas. While the results shown 
in Figure [4] emphasize the improvement in coherence, C- 
PACS also improves the visibility phases, which results 
in higher image fidelity. 

4.2. PP 13S* 

After applying the C-PACS atmospheric delay correc- 
tions to the 1.3-mm data, we removed the instrumental 
delay drifts by phase referencing to the 3C111 observa- 
tions. Figure p] shows the resulting maps before and after 
applying the (J-PACS corrections for the A and B con- 
figurations separately and the combined data sets. 

C-PACS improved the image quality for both the A 
and B configuration maps as measured by the increase 
in the peak flux, the reduction in the noise level, and 
the decrease in the observed source size. The improved 
image quality resulted from correcting the phase fluctua- 
tions in the 1.3 mm data. The increase in the source 
flux and decrease in source size is also illustrated in 
Figure |6j which shows radial profile plots across the 
1.3 mm emission toward PP 13S* along right ascension 
and declination. In the combined A+B configuration 
map, the peak flux measured toward PP 13S* increased 
from 42.4 mjy to 67.8 mJy (a factor of 1.6) after apply- 
ing the C-PACS correction. The noise level decreased 
from a — 1.5 mJy/beam to a = 1.1 mJy/beam, which 

12 The ^-distance is the baseline distance projected perpendic- 
ular to the line of sight. 



corresponds to a 36% improvement. The observed full- 
width-at-half-maximum (FWHM) source size diminished 
from 0.41" x 0.27" to 0.27" x 0.26", or a 52% decrease in 
the size of the major axis of the 1.3 mm emission. 

5. PROPERTIES OF THE PP 13S* CIRCUMSTELLAR DISK 

Before proceeding to the analyze the properties of the 
dust surrounding PP 13S*, we must flux calibrate the C- 
PACS data. The absolute flux calibration was set from 
observations of Uranus in B configuration and 3C84 in 
A configuration. The flux density of Uranus was inferred 
from a planet model, while the flux density of 3C84 was 
obtained from CARMA observations on a different day 
when both 3C84 and Uranus were observed. The un- 
certainty on this calibration is estimated to be 20% due 
to uncertainties in the planetary model and the boot- 
strapped flux for 3C84. The antenna gains as a function 
of time were then determined from the 3C111 observa- 
tions. 

Figure [7] shows the calibrated visibility amplitude ob- 
served toward PP 13S* as a function of baseline length. 
An unresolved source will have a constant flux density 
with baseline length. By contrast, the visibility ampli- 
tude toward PP 13S* decreases with increasing baseline 
length, which suggests that the source is resolved. While 
the decline in amplitude with baseline length could be 
explained by a coherence as low as ~ 0.2 on > 1-km 
baselines, the minimum measured coherence on 3C111 
at any ww-distance was 0.65 even before applying the 
C-PACS correction (see Figure H). A similar decline in 
the visibility amplitude for PP 13S* with baseline length 
can be seen as well when only a fraction of the data 
are averaged together (for example using 1 h of data at 
a time), indicating that atmospheric decorrelation over 
long timescales is not giving rise to the amplitude drop at 
long baselines. Thus, the primary cause of the decrease 
in amplitudes with increasing baseline length is that the 
source is resolved. 

The FWHM of the 1.3 mm continuum emission toward 
PP 13S* is 0.22" x 0.21", which was obtained by fit- 
ting a two-dimensional gaussian to the surface brightness 
distribution in the combined A+B configuration image 
after applying the C-PACS corrections and deconvolv- 
ing the synthesized FWHM beam size of 0.15" x 0.14". 
The integrated flux density of PP 13S* in the C-PACS 
corrected image after flux and amplitude calibration is 
241 ± 48 mJy, measured by integrating over an aperture 
of radius R dlsk ~ 1.4 x FWHM ~ 0.36", where R disk 
is defined to encompass 95% of the emission. The flux 
observed with CARMA corresponds to about half of the 
emission measured by single-dish obser yations (450 mJy 
at 1.3 mm with a beam size of 19.5" |Sandell fc Aspin 
19981. The remaining 1.3 mm flux is presumably con- 
tamed in an extended envelope larger than 1.6", which 
is the largest angular scale probed by the CARMA data. 

The presence of a circumstellar disk in PP 13S* has 
been previously inferred from several lines of evidence: 
(i) reflected light along the outflow axis is observed de- 
spite the fact that the centr al object is heavily obscured 
in the optical (Ay ~ 30—50; Cohen et al. 1983J), which in- 
dicates that the circumstellar material is not spherically 
symmetric; (ii) infrared absorption bands at 3/xm and 
10/im indicate substantial quantities of col d dust, prob- 
ably present in an obscured inclined disk (Cohen et al 
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1983| |Smith||1993[ ); and (iii) the broad 2.2 /an CO over- 
tone a bsorptionfe ature present in the PP 13S* spectra 



(Sandell & Aspin 1998 



Aspin & Sandell 2001) can be 



explained by thepresence ot a massive ac creting circum 



stellar disk (Hartmann & Kenyon 1996). Nonetheless, 
without observations ot the gas kinematics, we cannot 
determine if the dust emission detected by CARMA orig- 
inates from the central cusp of an envelope or from the 
circumstellar disk that surrounds PP 13S*. Since the 
presence of a massive accretion disk has been invoked to 
explain several characteristics of FU-Orionis objects, we 
assume that the millimeter continuum emission around 
PP 13S* observed by CARMA originates primarily from 
a circumstellar disk. 

To determine the disk properties, we assumed the ra- 
dial surface density [E(i?)] can be described by the sim- 
ilarity solution for a viscous accretion disk given by 



E(iJ,i) = S t 



exp 



1 



2(2- 7 ) 



2-7 



-(5) 



where 7 is the slope of the disk visc osity jv(R) oc IP) , 
and E t is the surface density at R t ( |Isella et al.| [2009). 
The transition radius is the radius at which the mass flow 
is zero, such that for R < R t the mass flow goes inward 
and mass is accreted into the disk, and for R > R t the 
mass flow goes outwards as the disk expands to conserve 
angular momentum. 

We assume that the ce ntral star has a bo lometric lu- 
minosity of Lb i = 30 L ( Cohen et al.|1983 ) and a mass 
of M* = 1 M Q . The dust opacity, assumed to be con- 
stant throughout the disk, is calculated using compact 
non- porous spheric a l grai ns with fractional abundances 



from |Pollack et al.| ( |1994[ ): 12% silicates, 27% carbona- 
ceous inaTeTIaTs^na 61% ices. The grain size distribution 
is assumed to be a power-law (n(a) oc a~ q ), with slope 
q = 3.5 and minimum grain size a m i n = 0.05 fxm. We 
adopt a dust emissivity index of /3 — 1.23 (derive d from a 
gray-body fit t o the spectral energy distribution; |Sandell| 
fc Aspin|1998[ ), and a dust-to-gas ratio of 0.01. For these 
assumed parameters, the implied maximum grain size of 
the distribution is a max = 0.1 cm and the mass opacity 
corresponds to k„ = 0.1 cm 2 g _1 at 1.3 mm. 

The dust surface density defined in Equation [5] was 
fitted to t he observed v i sibiliti es using the procedure de- 
scribed in |Isella et ai] ( |2009[ ). Figure [7] compares the 
observed and modeled visibility profile for PP 13S*. The 
model provides a reasonable fit to the data, although the 
observed visibility amplitudes are larger than the model 
for baselines longer than ~ 1 km. A map of the model 
was created and subtracted from the observed PP 13S* 
map, but no significant residuals (> 3tr) were found. 

The best fit disk model has a disk inclination of 15° 
(where 90° is defined as an edge-on disk), R t =13 AU, 
Et =145 g cm" 2 , 7 = 0.95, M disk =0.06 M , and 
i?disk=128 AU. The mass estimate is larger than the 



median mass of a Class II circumstellar disk (Andrews 
& Williams 2005), and is in agreement with the mea- 
sured masses around other FU-Orionis objects. For ex- 



ample Sandell & Weintraub (2001 ) estimated circumstel- 
lar masses (disk and envelope) between 0.02 M Q to a few 
solar masses for a sample of 16 FU-Orionis objects. 
Figure [8] shows the surface density distribution (E) and 



the optical depth (r) as a function of the disk radius. The 
vertical line marks the separation between the optically 
thin and thick regimes. We find that the disk becomes 
optically thick inwards of ~ 48 AU given the assump- 
tions in the model. We caution that the surface den- 
sity distribution within this region is poorly constrained 
not only because of the high optical depth, but because 
the disk is unresolved for a radius less than 26 AU. The 
radius where the dust becomes optically thick is large 
compared to what is found in disks around classical T 
Tauri stars, where only the inner few astronomical units 
are opaque at such long wavelength. Furthermore, the 
transition radius we find for PP 13S* is smaller than 
any of the pre-main sequ ence circumstellar disks stud- 
ied by llsella et al. (2009). Thus the circumstellar disk 
of PP 13S* is more concentrated than the disks around 
classical T- Tauri stars. 

Gravitational instabilities in the disk might be respon- 
sible for the en hanced accretion epis odes seen in FU- 
Orionis objects ( Armitage et al.pOOl ). To investigate if 
the disk around PP 13S* is gravitationally unstable, we 
computed the Toomre's Q parameter: 



Q = 



7rGE(i?) : 



(6) 



where c s is the sound speed and n is the epicyclic fre- 
quency (which is equal to the angular velocity for a kc- 
plerian disk). A disk becomes gravitationally unstable 
if Q 1-5, as spiral waves develop and mass is trans- 



Lodato & Rice 


2004) 


As shown m 1 


( 'igure 



§ the PP 13* disk is gravitation- 
ally stable to axisymmetric perturbations across all radii 
for the inferred surface density. The disk surface density 
would need to be increased by more than an order of 
magnitude (see dashed line of Figure [9| for the disk to 
develop a gravitational instability. The inferred disk sur- 
face density can vary widely depending on the adopted 
dust properties, as composition, grain size distribution 
(<Z, a m i n , a max ), emissivity index (/?) and dust-to-gas 
ratio, all affect the resulting opacity. For example, k v 
diminishes by 20-30% if ices are ignored from the dust 
composition, increasing the surface density by the same 
percentage. Also, flattening the grain distribution slope 
to q = 3.0 gives a 30% increase in the disk mass. Fur- 
thermore, we adopted a fixed power-law (/3) of the dust 
opacity law. If j3 decreases toward the center of the s ys- 
tem (as observed in Class sources; Kwon et al.|2009 | we 



can expect to have larger grains in the circumstellar disk 
that will reduce the mass opacity and increase the sur- 
face density. Despite the uncertainties in the dust prop- 
erties, it will be difficult to increase the surface density 
by more than an order of magnitude at a radius > 48 AU 
where the disk is optically thin, unless the dust proper- 
ties of PP 13S* are extr aordinarily different from what 
is found in typic al disks (Andrews & Williams 20 07a|b 
Isella et al.|[2009| 



6. CONCLUSIONS 

We have described C-PACS, which uses paired anten- 
nas as a means to calibrate the atmospheric phase fluc- 
tuations on long interferometric baselines. Specifically, 
while the 6 and 10-m CARMA antennas observe a sci- 
ence source in the 3-mm or 1-mm atmospheric windows, 
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the 3.5-m CARMA antennas simultaneously observe a 
nearby atmospheric calibrator in the 1-cm band. The 
3.5-m antennas are placed within 30 m of the larger an- 
tennas to sample similar atmospheric delay fluctuations. 
We have applied the calibration technique to CARMA 
observations of the circumstellar material around the 
FU Orionis object PP 13S*. C-PACS yields quantita- 
tive improvement in the image quality of PP 13S*: the 
observed peak flux increased by a factor of 1.6, the im- 
age noise level decreased by 36%, and the FWHM of the 
major axis decreased by 52%. 

The improvement in the phase error and amplitude co- 
herence provided by the paired antennas technique is a 
function of the projected beam separation of the antenna 
pairs at the height of the turbulent layer, and the radio- 
metric noise introduced from the reference array to the 
science array. Thus, the brightness of the atmospheric 
calibrator and the angular separation between the atmo- 
spheric calibrator and the science target are the main re- 
strictions to the application of this technique for general 
science observations. Our current estimate requires a cal- 
ibrator closer than 5° to the science target and brighter 
than 1 Jy at 30 GHz to correct 1.3 mm observations. 
Based on existing radio catalogs, we estimate that there 
are 420 sources that have S v > 1 Jy, such that 50% of 
the sky can be observed with a suitable calibrator. 

With C-PACS, we have obtained 0.15" resolution im- 
ages of the circumstellar material around PP 13S* at 
an observing frequency of 227 GHz. We measure an in- 
tegrated flux density of 241 mJy at 227 GHz, which is 



about half of the extended emission detected in a 19.5" 



beam ( Sandell k, Aspin||l998[ ). We constrain the surface 
density profile of PP 13S* using a self-consistent disk 
model. The main difference in the inferred disk proper- 
ties compared to disks around other pre-main sequence 
circumstellar disks is that the dust is more centrally con- 
centrated and there is a larger region that is optically 
thick at millimeter wavelengths. From analysis of the 
Toomre Q parameter, we find that the disk is gravita- 
tionally stable over all disk radii unless the disk surface 
density is underestimated by an order of magnitude or 
more. 
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Fig. 1. — Schematic of the CARMA Paired Antenna Calibration System. Antennas Ai and A2 (in the science array) observe the science 
target while antennas A3 and A4 (in the reference array), offset by a distance b from the science array, observe a bright atmospheric 
calibrator. In the turbulent layer at an elevation h the beams are separated by a distance b e ff that depends on b, h, and the angular offset 
between the science and reference sources, a. If the effective baseline, b B g, is much smaller than the science baseline, B, the delay difference 
to the reference antennas is a good estimate for the delay difference to the science antennas. 
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Fig. 2. — The left panel shows the locations of the science (10.1-m and 6.4-m) and reference (3.5-m) antennas in the B-configuration 
of CARMA (circles do not represent the actual antenna diameters). The right panel shows the separation, b e g, between the science and 
reference beams in the turbulent layer for each pair of antennas during the PP 13S* observations. The continuous (dashed) line corresponds 
to a turbulent layer at a l(2)-km altitude. For most pairs the configuration is very favorable, giving rise to effective baselines much shorter 
than the science baseline. 
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Fig. 3. — Visibility phase versus time for the science array (orange) and the reference array (green) measured toward 3C111 on a baseline 
of length ~ 800 m in B configuration. The phase measured by the reference array at a frequency of 30.4 GHz was scaled to the frequency 
of the science array (227 GHz) by a factor of 7.5 (= 227 GHz / 30.4 GHz). Each data point indicates the measured phase over a 4 s 
integration after removing the mean phase computed in 10 min intervals. About 1 h of data are omitted near the middle of the observation 
when PP 13S* and 3C111 were transiting at an elevation > 80° and the tracking of the antennas was poor. The figure demonstrates that 
the observed phases at 30.4 GHz closely track the phases measured at 227 GHz and can be used to correct the atmospheric fluctuations at 
higher frequencies. 
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Fig. 4. — Visibility phases versus Jiti-distance on all paired baselines for observations of 3C111 in the CARMA B and A configurations. 
The phases are shown before (top panels) and after (middle panels) applying C-PACS corrections. Each point indicates the measured phase 
in a 4 s integration after removing the mean phase computed in 10 min intervals. The bottom panels show the coherence calculated over 
130 kA intervals. Before applying C-PACS, the coherence declines with iiu-distance as expected for atmospheric phase fluctuations that 
increase with baseline length. After applying C-PACS, the phase coherence is higher and uniform with baseline length. 



10 



Perez et al. 



Without C-PACS 



With C-PACS 



B configuration 



o 

<D 
en 
O 

'JO 
< 



o 

0) 
o 

< 



o 

w 
o 

< 




o-1.7 mJy/beam 



0-1.5 mJy/beam 




0.4 0.2 -0.2-0.4 

Ao [arcsec] 



0.4 0.2 -0.2-0.4 

Aa [arcsec] 



Fig. 5. — CARMA 227 GHz continuum images of PP 13S* before (left) and after (right) C-PACS corrections for data obtained in B 
configuration (top), A configuration (middle), and B and A configurations combined (bottom). No absolute flux scale or amplitude gain 
calibration have been applied to the data in order to assess the effect of atmospheric phase corrections only. The color scale range is same 
for all maps (from -2.4 to 54.8 mjy) such that the measured fluxes can be compared directly. Solid contours are at 2<r, 5a and in increments 
of 5a thereafter. Dotted contour is at —2a. For both B and A configurations, applying the C-PACS correction increased the observed peak 
flux and reduced the observed source size. 
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Fig. 6. — Observed flux density toward PP 13S* versus offset in right ascension (left) and declination (right) before and after applying 
C-PACS corrections. Applying C-PACS corrections increased the peak flux density and decrease the observed source size. 
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Fig. 7. — Observed visibility amplitudes toward PP 13S* after applying the C-PACS corrections (black points) versus the deprojected 
baseline length. The errorbars indicate the la interval uncertainties, but exclude uncertainties in the flux calibration (~20%). The 
histogram shows the expected signal for random noise. The dashed line shows the best-fit disk model (see text) to the visibility data. The 
clear decline in the visibility amplitude with baseline length indicates the dust emission around PP 13S* has been resolved. 
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Fig. 8. — Model surface density distribution of PP 13S* as a function of disk radius. The disk optical depth (r = Sd(-R) X rei.3 mm ) was 
computed assuming a constant dust opacity throughout the disk. The inner region of the disk (R < 48 AU) becomes optically thick at a 
frequency of 227 GHz. 
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Fig. 9. — Toomre's instability parameter, Q, plotted for the PP 13S* disk. The solid line shows the Q values for the nominal disk model. 
Dust opacities that are 10, 20 and 40 times larger than the adopted opacities are shown by the dashed, dash-dotted, and vertical-dashed 
curves respectively. The Q = 1.5 limit is shown for reference, such that regions with Q < 1.5 may become gravitational unstable. The disk 
surface density needs to be an order of magnitude higher than the nominal model in order to form instabilities in the disk. 



